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ABSTRACT 

We identified and characterized the largest (536) RR Lyrae (RRL) sample in a 
Milky Way dSph satellite (Sculptor) based on optical photometry data collected over 
~24 years. 

The RRLs display a spread in R-magnitude (~0.35 mag) which appears larger 
than photometric errors and the horizontal branch (HB) luminosity evolution of a 
mono-metallic population. Using several calibrations of two different reddening free 
and metal independent Period-Wesenheit relations we provide a new distance estimate 
^=19.62 mag ((7^=0.04 mag) that agrees well with literature estimates. We constrained 
the metallicity distribution of the old population, using the Mj Period-Luminosity 
relation, and we found that it ranges from -2.3 to -1.5 dex. The current estimate is 
narrower than suggested by low and intermediate spectroscopy of RGBs (A [Fe/H]< 
1.5). 

We also investigated the HB morphology as a function of the galactocentric dis¬ 
tance. The HB in the innermost regions is dominated by red HB stars and by RRLs, 
consistent with a more metal-rich population, while in the outermost regions it is 
dominated by blue HB stars and RRLs typical of a metal-poor population. Our re¬ 
sults suggest that fast chemical evolution occurred in Sculptor, and that the radial 
gradients were in place at an early epoch. 
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1 INTRODUCTION 

Understanding the evolution of Local Group dwarf galax¬ 
ies offers a fundamental key to constrain the early forma¬ 
tion and evolution of cosmic structures (e.g. Mayer 2010; 
Madau & Dickinson 2014). Current ground-based and space 
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facilities allow us to resolve their stellar content and to pro¬ 
vide a wealth of observables to reconstruct their star forma¬ 
tion history and chemical enrichment (Gallart et al. 2005; 
Tolstoy et al. 2009). 

RR Lyrae (RRL) stars are reliable distance indica¬ 
tors, since their individual distances can be estimated us¬ 
ing either the different developments of the Baade-Wesselink 
method (Storm et al. 2004), the correlation between visual 
magnitude and Fe abundance (Cacciari & Clementini 2003), 
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statistical parallaxes (Dambis et al. 2014), optical, near- 
infrared (NIR) and mid-infrared (MIR) Period-Luminosity- 
Metallicity (PLZ) (Bono et al. 2001, 2003; Catelan et al. 
2004) and the reddening free Period-Wesenheit-Metallicity 
(PWZ, Braga et al. 2015; Coppola et al. 2015). The above 
diagnostics all have pros and cons, but during the last few 
years the use of accurate and precise NIR and MIR mean 
magnitudes has provided a new spin on the use of RRLs as 
distance indicators (Madore et al. 2013; Klein et al. 2014; 
Neeley et al. 2015). This evidence was soundly comple¬ 
mented by the recent estimates of the trigonometric paral¬ 
laxes for hve field RRLs using the FGS on board of HST 
(Benedict et al. 2011). Moreover, the use of the redden¬ 
ing free magnitudes, the so-called Period-Wesenheit relation 
(PWR) is hrmly supported by theory (Marconi et al. 2015) 
and observations (Braga et al. 2015). 

RRL stars also play a fundamental role as stellar trac¬ 
ers. Theory and observations indicate that they are present 
in all stellar systems hosting an old (> 10 Gyr) stellar pop¬ 
ulation. They are low-mass (0.6 — 0.8M©) central helium¬ 
burning stars in their Horizontal Branch (HB) evolution¬ 
ary phase. This ancient population is the fossil record of 
the early stages of galaxy evolution and provides firm con¬ 
straints on the time-scale of their early formation and evo¬ 
lution. 

In this context. Sculptor is an interesting labora¬ 
tory, since it hosts a conspicuous old stellar population, 
but its star formation is quite complex (Da Gosta 1984; 
Majewski et al. 1999; Tolstoy et al. 2004; de Boer et al. 
2012). Moreover, Sculptor appears to be characterized by a 
complicated chemical enrichment history (Smith & Dopita 
1983; Majewski et al. 1999; Hurley-Keller et al. 1999; 
Starkenburg et al. 2013). In particular, the presence of a 
large metallicity spread has been suggested in this galaxy. 
Detailed spectroscopic measurements (Tolstoy et al. 2004; 
Walker et al. 2007, 2009) showed not only that the more 
metal-rich red giant branch (RGB) stars ([Fe/H]> -1.7) 
are more centrally concentrated than more metal-poor ones, 
leading to a significant metallicity gradient, but also that 
they have different kinematics. More importantly, the fact 
that red HB (RHB) stars are more centrally concentrated 
than blue HB (BHB) stars suggests that a spread in metal¬ 
licity is also present in the old stellar component (HB stars), 
which in turn suggests a relatively fast early chemical evo¬ 
lution. Since the stars that are now RR Lyrae stars were 
formed during an even shorter time than the global HB 
population, observing a gradient in their pulsational proper¬ 
ties would put even more stringent constraints on the time- 
scale of the chemical enrichment in Sculptor (Bernard et al. 
2008). 

Since the discovery of Sculptor by Shapley (1938), a 
number of estimates of its distance have been proposed in 
the literature. After the initial approximate value by Shapley 
(1938) based on photographic data, the first accurate mea¬ 
surement dates back to the work of Hodge (1965) who de¬ 
rived a true distance modulus of 19.7±0.15 mag using three 
different methods (luminosity of the RGB, two W Virginis, 
and three RRL stars). Kunkel & Demers (1977) provided 
distance estimates based on the average magnitude of 24 
stars on the HB with colours covering the RR Lyrae gap 
(0.15<(B-F)o <0.52); they found a true distance modulus 
of 19.47±0.10 mag, slightly lower than the previous one. 


More recently, Pietrzynski et al. (2008) provided a dis¬ 
tance value of /i=19.67±0.14 mag by applying different the¬ 
oretical and empirical calibrations of the PLZ relation for 
RRLs in the NIR (J, K) magnitudes, consistent with the 
F-band data of RRLs and also in very good agreement with 
the results obtained by Rizzi (2002) based on the tip of RGB 
(TRGB, 19.64T0.08 mag) and HB stars (19.66±0.15) based 
on optical photometry. 

The search for variable stars in Sculptor dates back 
to more than half a century ago (Baade & Hubble 1939; 
Thackeray 1950). The first systematic search was performed 
by van Agt (1978), who identified more than 600 candidate 
variable stars, but provided periods only for ~10% of the 
sample. 

The most complete RRL catalogue to date was provided 
by Kaluzny et al. (1995) based on the OGLE-I survey ob¬ 
servations (Udalski et al. 1992). They investigated the cen¬ 
tral region of Sculptor (15'xl5') and identified and char¬ 
acterized 226 RRLs. Their pulsation properties are consis¬ 
tent with metal-poor chemical composition ([Fe/H]<-“1.7) 
and a spread in metallicity. The spectroscopic follow-up by 
Glementini et al. (2005) confirmed this preliminary evidence 
based on low resolution (LR) spectra of 107 RRLs. In par¬ 
ticular, they found that the metallicity distribution peaks at 
[Fe/H]~ -1.8, but covers more than 1.5 dex (-2.4<[Fe/H]<- 
0 . 8 ). 

We take advantage of a large sample of RRLs, including 
some new discoveries, to provide a revised distance estimate 
for Sculptor and to constrain the parameters of an early 
chemical enrichment. The current sample is a factor of two 
larger than previously analysed ones (536 vs 216^). More¬ 
over, we provide pulsation period, mode and mean magni¬ 
tudes in three different photometric bands (see § 2). In § 4 we 
describe the approach for the new distance determination, 
while in § 5 we investigate the spread in metallicity of the 
old stellar population using the /-band Period-Metallicity 
relation of RRLs. 


2 PHOTOMETRIC DATA SET AND RR 
LYRAE IDENTIFICATION 

We used 4,404 calibrated images from the photometric 
database collected by one of us (PBS), covering an ex¬ 
tended area around the Sculptor dSph. The time base¬ 
line of these data is nearly 24 years (from October 1987 
until August 2011). These data are individual CCD im¬ 
ages from many ground-based telescopes: CTIO 0.9 m-|-[TI, 
Tek2K], LCO 1.0 m-bFORD2 (OGLE-I survey), MPI/ESO 
2.2m-bWFI, ESO NTT 3.6m-bSUSI, CTIO 4m-b Mosaic2, 
SOAR 4.1 m-|-SOI. The full data-set covers an area of -^6 
sq. deg. centred on the Sculptor dSph galaxy. However, 
only the central region (-^3 sq. deg) could be properly cali¬ 
brated. The DAOPHOT/ALLFRAME package of programs 
(Stetson 1987, 1994) was used to obtain the instrumental 
photometry of the stars. A detailed description of the obser¬ 
vations and data reduction will be given in a forthcoming 
paper. 

^ 10 out of the 226 RRLs published by Kaluzny et al. (1995) 
turned out to be either double identifications or non-variables 
stars. 
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We performed the variability search over the full data 
set. An updated version of the Welch-Stetson variability in¬ 
dex (Welch & Stetson 1993) was used to identify candidate 
variable stars on the basis of our multi-band photometry. 
From the list of 663 variable star candidates, we have found 
611 actual variable stars and 23 possible variable stars. Out 
of these real variables, 536 pulsating stars are located in the 
Instability Strip (IS) of a very extended HB of the Sculptor 
dSph galaxy, with periods between 0.2 and 1.2 days. Out of 
these, 82 RRLs are new identifications, and for 320 RRLs all 
their pulsational parameters (period, mean magnitude, and 
amplitude) are derived for the first time. This work increases 
by over a factor of two the number of the known RRL to 
date in the Sculptor dSph. Our sample includes 216 RRLs 
of Kaluzny et al. (1995). 

We derived pulsational properties for all the RRL stars 
from our R F/-Johnson/Cousins photometry. The search for 
the period was carried out using a simple string-length al¬ 
gorithm (Stetson et al. 1998). Then a robust least-squares 
fit of a Fourier series to the data refined the periods. The 
intensity-averaged magnitudes and amplitudes of the mono- 
periodic light curves were obtained by fitting the light 
curves with a set of templates partly based on the set of 
Layden et al. (1999) following the same method described 
in Bernard et al. (2009). Thus, through the period and 
light-curve shape, we identified: i) 289 RRab, pulsating in 
the fundamental mode; 20 of which are suspected Blazhko, 
(Blazhko 1907)); ii) 197 RRc, pulsating in the first-overtone 
mode; and in) 50 possible multi-mode RRd stars, pulsating 
in both modes simultaneously, although the classification in 
some cases was uncertain due to their relatively noisy or 
(very) poor light curves. The mean (maximum) number of 
points in the light-curves of the RRL stars are 83, 52, and 
21 (115, 182 and 28) respectively in B, V, and I. 


3 RR LYRAE INSTABILITY STRIP 

The top left panel Fig. 1 shows the (B~V, V) colour- 
magnitude diagram (CMD) of Sculptor zoomed on the HB 
region. Data plotted in this CMD show a complex mor¬ 
phology with well defined BHB and RHB stars (black dots) 
and a large sample of RRLs. The latter group shows sev¬ 
eral interesting features, i) The RRL IS is well populated 
over the entire colour range typical of RRc and RRab vari¬ 
ables.the spread in magnitude of RRLs is ~0.35 mag, 
significantly larger than the typical uncertainty in the mean 
magnitudes ((7=0.03 mag) and significantly larger than the 
expected from the evolution of a mono-metallic population. 
The top-right panel of Fig. 1 shows the visual magnitude 
distribution of the RRLs located within 2.5 a of the aver¬ 
age. We ended up with a sample of 520 RRLs and neglected 
16 outliers. Using the period-amplitude distribution and by 
visual inspection of the individual light curves we identified 
276 RRab + 195 RRc (defined by us as full sample) and 49 
candidate RRd variables, iii) The magnitude distribution of 
the RRLs is broad and possibly double peaked (see top-left 
panel). 

A plausible working hypothesis to account for the 
observed spread in magnitude of RRLs is that Sculptor 
RRLs cover a broad range in metallicity. Using classical 
luminosity-metallicity relations My {RRL) = a-\- j3 ■ [Fe/H] 


by Clementini et al. 2003 and Bono et al. 2003, we found 
that a spread in U-magnitude of 0.35 mag implies a range 
in Fe abundance of the order of ~0.6 dex. Note that this 
range has to be treated with caution, since the above rela¬ 
tions are prone to systematic uncertainties, and in particular 
to evolutionary effects (Bono et al. 2001, 2003). 

In order to search for possible metallicity trends, we 
split the sample of RRab and RRc stars into two groups - 
candidate RRd might have less accurate mean magnitudes 
and were thus neglected -. We have used the mean mag¬ 
nitude of the entire sample, <U>=20.13 ((7=0.09) mag, as 
the limit, shown by a solid line on the right-panel and the 
arrow on the right-panel of Fig. 1. We note that this selec¬ 
tion is marginally affected by the curvature of the HB, in 
particular for the bluer stars, but that the results did not 
change significantly with a different criterion. We find 184 
bright (Bt; blue symbols) and 287 faint (Ft; orange symbols) 
RRLs. 

A more restrictive selection based on the quality of the 
phase coverage of the photometry over the entire pulsation 
cycle yielded similar results. This refined sample contains 
290 RRLs that we defined as the clean sample. Among them 
we have 167 RRab and 123 RRc variables. Data plotted in 
the bottom panels of Fig. 1 indicate that the global prop¬ 
erties are quite similar to the full sample, and indeed the 
mean V magnitude is < 1/>=20.14((7=0.07) mag. The same 
conclusion applies to the mean V magnitude of both Bt- 
and Ft-RRLs (see labelled values). 

In order to rule out the hypothesis that evolutionary 
effects can explain the observed luminosity range. Figure 
2 shows a comparison with stellar evolutionary tracks for 
different metallicities. The current predictions come from 
the BaSTI library (Pietrinferni et al. 2004), they have been 
computed assuming a scaled-solar chemical mixture (using 
the cr-enhanced HB models provide consistent results). 
They were plotted assuming the reddening law from 
Cardelli et al. (1989) and assuming a true distance modulus 
of /i = 19.62 mag (see § 4) and a reddening of 0.018 mag. 
The three panels present the case for Z=0.0001, 0.0003, 
0.0006. In each panel the red thick line shows the Zero Age 
Horizontal Branch (ZAHB), while individual evolutionary 
tracks for selected labelled masses are presented by the thin 
black lines. The black circles mark the points of each tracks 
with a time step of 10 Myr after the He ignition. Overall, 
the figure shows many interesting points: 

i) The luminosity of HB stars mainly depends on the 
metallicity, therefore it is plausible to assume that the faint 
sample previously selected includes the most metal-rich 
RRL stars, while the bright sample should include the more 
metal-poor stars. However, off-ZAHB evolution could move 
more metal-rich stars into the bright sample. 

ii) The comparison with tracks with Z=0.0006, in the right 
panel, suggests that this is a fair upper limit to the RRL 
metallicity. Nevertheless, the ZAHB is not able to well 
reproduce the RHB, suggesting that even more metal-rich 
stars populate this feature of the CMD. 

iii) The left panel (Z=0.0001) shows that very metal-poor 
stars are expected to start core Helium-burning on the 
blue side of the IS, except possibly the most massive ones. 
Nevertheless, the stars evolving off the ZAHB do cross 
the IS, although this occurs at relatively bright luminosity 
where only a few very bright RRL stars are observed. 
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Figure 1. Top-left. Optical {B-V, V) CMD of Sculptor, zoom-in on HB stars. Blue and orange symbols display bright and taint RRLs 
of the full sample. The horizontal line shows the mean magnitude of the full sample (1/=20.13 mag) adopted to split the bright and 
faint samples. Diamonds and circles display RRab and RRc variables. The number of RRLs in each subsample is labelled. Bottom-left. 
Same as the top, but for the clean sample. Top-right. F-band luminosity distribution of RRLs. The arrow marks the magnitude adopted 
to split bright and faint RRLs. Bottom-right Same as the top, but for the clean sample. For the sake of clarity, RRd are not represented. 


Moreover, the spacing between black points indicates 
that the crossing of the IS is an extremely rapid phase, 
thus very difficult to observe. This is expected to be even 
more effective for extremely metal poor stars (Z<0.0001). 
Therefore, although possible, it seems unlikely that many 
of the Sculptor RRL stars are as metal-poor as the most 
metal-poor stars detected in the RGB (Starkenburg et al. 
2013). 

iv) While the tracks corresponding to each metallicity 


cover a significant range of luminosity, it is clear that the 
luminosity spread associated to the evolutionary effects 
of a mono-metallic population is not sufficient to explain 
the observed magnitude range of RRL stars. This strongly 
supports our hypothesis that the RRLs in Sculptor cover a 
wide range of metallicities. 
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Figure 2. Optical (V,B-V) CMDs of Sculptor, zoom-in on HB stars. The grey circles display the RRL variable stars of the full sample 
(same as top-left panel in Fig. 1). The red lines show the zero age horizontal branch (ZAHB) from the BaSTI library (Pietrinferni et al. 
2004) for different metallicities: Z=0.0001 (left), Z=0.0003 (middle), and Z=0.0006 (right). The black lines display the corresponding 
evolutionary HB models with the labelled masses, where the overplotted circles indicate time intervals of 10 Myr. The stellar mass, at 
the typical colours of the R,R, Lyrae stars, steadily increases when moving from the brightest to the faintest HB models. 


Table 1. Parameters of the theoretical fit, W(V,B-I) = a + fi 
■ logP[mag], performed following the procedure in Marconi et al. 
(2015) for metal-independent Wesenheit magnitudes. 


TH W(V,B-I) 

a. 

p 

cr 

RRc 

-1.78±0.06 

-3.12±0.13 

0.06 

RRab 

-1.134±0.008 

-2.60±0.04 

0.07 

Global 

-1.134±0.009 

-2.48±0.04 

0.08 


4 DISTANCE DETERMINATION 

The evidence that the RRLs in Sculptor show a spread 
in chemical composition is a thorny problem for distance 
determinations. The diagnostics adopted to estimate indi¬ 
vidual RRL distances depend on the Fe content (see § 1). 
However, recent theoretical (Marconi et al. 2015) and em¬ 
pirical (Braga et al. 2015; Coppola et al. 2015) evidence in¬ 
dicates that the Period-Wesenheit relation (PWR) in the 


{V, B-V) and the {V, B-I) bands are minimally affected 
by metallicity. The above Wesenheit magnitudes are de¬ 
fined as W{V, B-V)=V-3.06{B-V) and W{V, B-I)=V- 
1.34(B-/), where the coefficients were fixed according to the 
extinction law provided by Cardelli et al. (1989) calculated 
for the effective mean wavelengths of Landolt’s (Landolt 
1992) filter/photocathode combination for stars of spectral 
type A0III-F5III. 

In their investigation Marconi et al. (2015) provided a 
detailed analysis of optical, optical-NIR and NIR Period- 
Wesenheit-Metallicity (PWZ) relations of the form W{X, 
X“Y) = a + P ■ logP-h 7 • [Fe/H], where X and Y are mag¬ 
nitudes. They found that the coefficient of the metallicity 
term was smaller than ~ 0.05 dex for the W{V,B-V) and 
the W{ V, B-I) PWR (see their Tables 7 and 8), indicating 
that the metallicity dependence of the F-band is counter¬ 
balanced by the convolution between the reddening vector 
and {B-V)-{B-I) colors. 


MNRAS 000, 1-10 (2015) 













6 C. E. Martinez-Vazquez et al. 



- 0.5 - 0.4 - 0.3 - 0.2 - 0.1 - 0.5 - 0.4 - 0.3 - 0.2 - 0.1 

log P log P 


Figure 3. Top-left. Period-Wesenheit {V,B—V) relation for RRab and RRc RRLs of the clean sample. Symbols are the same as Fig. 1. 
Black lines display the least squares fits of the two subsamples. The sigma of the fits are labelled. Top-right. Same as the left, but for 
the global {RRab-\-RRc) sample fundamentalized. Bottom-left. Same as the top left, but for the PW ( V, B-I) relation. Bottom-right. 
Same as the bottom left, but for the global RRL sample fundamentalized. 


The top-left panel of Fig. 3 shows the PWR( V, B- V) of 
the clean RRL sample. As expected RRc and RRab variables 
display different PWRs. The vanishing dependence of the 
quoted PWRs on metallicity is further supported by the 
small dispersion in apparent Wesenheit magnitude of RRLs 
for each period. We performed a linear least squares fit to the 
three different subsamples and also applied a S-a clip, and 
the final linear regressions are plotted as black lines in the 
top panels of Fig. 3. The coefficients of the empirical PWRs 
are listed in Table 2. To constrain the impact of using the 


clean sample on the empirical PWRs, the least squares fits 
were also performed using the full RRL sample. Data listed 
in Table 2 clearly indicate that the coefhcients of the PWRs 
attain similar values within the uncertainties. 

We can see how the Bt-RRLs and Ft-RRLs follow the 
same relation. Moreover, Bt-RRLs seem to be the natural 
extension in the long-period regime of Ft-RRLs. The large 
sample of RRLs variables allows us to estimate the distance 
to Sculptor using three different PWRs for RRc, RRab, and 
the global sample {RRcf-RRab, right panel of Fig. 3). In 
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Table 2. Parameters of the empirical fit, W= a + /3 ■ logP[mag], performed on each group of RRL stars. 


EM 


CLEAN 



FULL 


RRc 

a 

p 

<T 

OL 

p 

CT 

W{V,B-V) 

W{V,B-I) 

18.21±0.06 

17.95±0.05 

-2.61±0.13 

-2.92±0.10 

0.06 

0.04 

18.26±0.05 

18.05±0.04 

-2.48±0.10 

-2.65±0.08 

0.06 

0.05 

RRab 

a 

p 

<7 

OL 

p 

cr 

W{V,B-V) 

wlv,B-I) 

18.56±0.03 

18.38±0.02 

-2.66±0.12 

-2.88±0.07 

0.08 

0.05 

18.54±0.02 

18.38±0.02 

-2.72±0.10 

-2.83±0.07 

0.09 

0.06 

Global 

a 

p 

<7 

OL 

p 

CT 

W(V,B-V) 

W(V,B-I) 

18.59±0.02 

17.42±0.01 

-2.48±0.05 

-2.66±0.04 

0.07 

0.05 

18.58±0.01 

18.44±0.01 

-2.50±0.05 

-2.54±0.03 

0.08 

0.06 


Table 3. Distance moduli (in mag.) based on the current sample of RRL. 





CLEAN 



FULL 


RRab 


M (V, B-V) [165] 

Ai (U, B-I) [164] 

< fj. > 

Al (V, B-V) [261] 

Al (y, B-I) [257] 

< fj. > 

TH 

SE 

EM (o = 

[-1.03±0.10; -1.04±0.10]) 

19.64±0.02(0.08) 

19.64±0.02(0.08) 

19.59±0.02(0.08) 

19.57±0.01(0.05) 

19.51±0.01(0.05) 

19.59(0.05) 

19.63±0.02(0.09) 

19.61±0.02(0.09) 

19.58±0.02(0.09) 

19.57±0.01(0.06) 

19.52±0.01(0.06) 

19.58(0.04) 

RRc 


M (V, B-V) [116] 

Ai (U, B-I) [114] 

< fj. > 

Al (V, B-V) [178] 

M (V, B-I) [177] 

< fj. > 

TH 

SE 

EM (o = 

[-1.31±0.17; -1.24±0.17]) 

19.70±0.02(0.06) 

19.77±0.02(0.06) 

19.52±0.02(0.06) 

19.63±0.01(0.04) 

19.73±0.01(0.04) 

19.67(0.10) 

19.69±0.02(0.06) 

19.82±0.02(0.08) 

19.51±0.03(0.06) 

19.62±0.01(0.05) 

19.84±0.01(0.05) 

19.70(0.14) 

Global 


M (V, B-V) [279] 

Ai {V, B-I) [277] 

< P > 

Al (V, B-V) [437] 

Al (V, B-I) [429] 

< p > 

TH 

SE 

EM (o = 

[-0.98±0.09; -0.98±0.09]) 

19.66±0.02(0.07) 

19.66±0.02(0.07) 

19.57±0.02(0.07) 

19.60±0.01(0.05) 

19.55±0.01(0.05) 

19.61(0.05) 

19.65±0.02(0.08) 

19.65±0.02(0.08) 

19.56±0.02(0.08) 

19.60±0.01(0.06) 

19.58±0.01(0.06) 

19.61(0.04) 

< M > ADOPTED 



19.62(0.04) 



19.63(0.06) 


deriving the distances of the global sample the periods of the 
RRc variables are fundamentalized, i.e., logRp = log Pfo + 
0.127 (Bono et al. 2001; Marconi et al. 2015). 

The bottom panels of Fig. 3 display the PWR for RRLs 
in Sculptor, but in the V, B-I bands (see Table 2). The 
standard deviations of the three different PWRs are slightly 
smaller than in V, B- V. The difference is mainly due to 
the fact that B~I is a better temperature indicator and 
the PWR closely mimics a Period-Luminosity-colour rela¬ 
tion (Bono et al. in preparation). 

The distance estimates to Sculptor were derived using 
the three different subsamples {RRab, RRc, and global) and 
both the clean, 290 and the full, fll RRL stars. In order to 
quantify the effect of both the zero-point (a) and the slope 
(/3) of the relation on the distances, we adopted three dif¬ 
ferent approaches: i) theoretical - TH- in which we adopted 
predicted PWRs (Marconi et al. 2015); ii) semi-empirical - 
SE- in which we adopted the observed slopes listed in Ta¬ 
ble 1 and the theoretical zero-points; Hi) empirical - EM - in 
which we adopted the empirical slopes and where the zero- 
points are based on the field RRLs with HST trigonometric 
parallaxes (Benedict et al. 2011). In particular, we adopted 
RR Lyr itself to calibrate the RRab and the global sample, 
while the RRc were calibrated using RZ Cep. The mean 
B and V magnitudes provided by Coppola et al. (2015) 
and the individual distance moduli provided by Braga et al. 
(2015) and Neeley et al. 2015 have been used. 

Using the empirical slopes and the above distances we 


estimated the zero-points for RRab, global and RRc PW re¬ 
lations (see Table 3). The former values for RRab and global 
relations agree quite well with predicted ones (-1.08±0.01, 
-1.07±0.01). However, in the case of the RRc, the new zero- 
points are ~30% smaller than predicted (-1.56±0.04). 

Table 3 gives the distance moduli based on the clean 
and full RRL samples interpreted with the three different 
PWRs {RRab, RRc, global) and the three different calibra¬ 
tions (TH, SE, EM). The number of RRLs used to estimate 
the distance after the 3a clipping is given in square brack¬ 
ets. Table 3 gives true distance moduli (mean values of p 
( V, B-V) and p (V, B-I) over the entire sample. For each 
fixed sample and approach, the estimated distance moduli 
agree with each other within Icr and further support the ev¬ 
idence that the two clean and full samples give very similar 
distances. Since there is no reason a priori to prefer one of 
the five different estimates, for each sample we estimated the 
mean value —denoted hy < p >— that is given with the 
relative standard deviation (cr^). The distance determina¬ 
tions listed in Table 3 agree quite well with similar distances 
available in the literature. In particular, they agree with the 
distance obtained by Pietrzynski et al. (2008) using RRLs 
[19.67±0.02(0.12)]. The difference is indeed smaller than 
0.4a. The true distance modulus to Sculptor was estimated 
as the mean among the mean distances of RRc, RRab and 
global sample. The mean distance based on the global sam¬ 
ple is treated as an independent determination, because the 
RRc, once fundamentalized, become “pseudo” fundamental 
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Figure 4. Top. Absolute Mj Period-Luminosity relation for 
RRab and RRc variables of the clean sample. Symbols are the 
same as in Fig. 1. The dashed and solid black lines show pre¬ 
dicted PLR for two different Fe abundances ([Fe/H]=-1.5, -2.3; 
(Marconi et al. 2015)). Bottom. Same as the top, but for the 
global sample. 

pulsators (Braga et al. 2015). We found < fi >cZean=19.62, 
cr^=0.04 and < /r >/uii=19.63, cr^^O.OB mag, respectively. 
We adopted the former one, for the reasons mentioned 
above. 

5 METALLICITY OF THE RRL STARS 
5.1 Metallicity range 

We will now take advantage of the new accurate distance de¬ 
termination to Sculptor [19.62±0.01(0.07)] and the depen¬ 
dence of the /-band Period-Luminosity (PLI) relation on 
metallicity (Marconi et al. 2015) to constrain the metallic¬ 
ity distribution of the clean RRL sample. The top-left panel 
of Fig. 4 shows the comparison between RRLs and predicted 
PLI relations for two different metal abundances [Fe/H]=- 
2.3 (solid line) and [Fe/H]=-1.5 (dashed line). The majority 
of both RRc and RRab are bracketed by the quoted pre¬ 
dictions. This evidence is further supported by the global 
sample (bottom-left panel). Moreover, the Bt-RRL and Ft- 
RRL subsample are located closer to the metal-poor and 
metal-rich relations, respectively. We performed the same 
test using the full sample and found that the metallicity 
ranges from ~ -2.4 to ~ -1.4 dex. 

The spread in metallicity of the RRLs that we have 
measured using this method - 0.8 to 1.0 dex - is in excellent 
agreement with that obtained by Clementini et al. (2005) 
from spectroscopic measurements of RRL stars, but signif¬ 
icantly smaller than the estimates based on spectroscopic 
metallicities of RGB stars. Measurements from LR spectra 


of 107 RRLs (AS parameter) indicate a mean Fe abundance 
of-1.83±0.03 (a—0.26) with the metallicity ranging from ~ 
-2.40 to -0.85 dex (Clementini et al. 2005). However, they 
found only one metal-rich ([Fe/HJ^-^ -0.85) RRL star while 
the bulk covers in a metallicity range from ^ -2.4 to - 
1.4 dex. Therefore, our estimates are in agreement with the 
bulk of the RRLs included in Clementini et al. (2005). More 
detailed spectroscopic measurements from medium resolu¬ 
tion spectra for more than 300 RGBs revealed the occur¬ 
rence of a metal-poor stellar component with Fe abundance 
ranging from -2.8 to -1.7 and of a metal-rich component 
with metallicity ranging from -1.7 to -0.9 (Tolstoy et al. 
2004; Coleman et al. 2005). The empirical scenario was re¬ 
cently fortified by the extensive spectroscopic investigation 
based on medium resolution spectra of almost 400 RGBs by 
Kirby et al. (2009). They found a mean metallicity of-1.58 
((t= 0.41) with the Fe content ranging from ^ -3 to ~ -1 
dex. Moreover, they also found evidence of a linear metal¬ 
licity gradient. 

However, we cannot exclude the possible occurrence of 
small samples of more metal-poor or metal-rich old stars, 
because these stellar components might produce HB stars 
either too hot or too cool to inhabit the RRL IS. This par¬ 
ticularly means that RRLs are not sensitive to the pos¬ 
sible occurrence of an extremely metal-poor component 
(Kirby et al. 2009; Starkenburg et al. 2013; Simon et al. 
2015, see also Fig. 2). However, the lack of High Amplitude 
Short Period fundamental RRLs (HASP; Log P < -0.35, 
Av > 1, Fiorentino et al. 2015) suggests that Sculptor does 
not host a significant metal-rich ([Fe/H]>-1.3) old stellar 
population (Martinez-Vazquez et al. in preparation). 

5.2 Radial properties 

Tolstoy et al. (2004) found a radial gradient in the metallic¬ 
ity distribution of both RGB and HB stars, in the sense that 
the more metal-poor component is more spatially extended 
than the metal-rich one. Now, we will explore whether a ra¬ 
dial gradient exists in the properties of the Sculptor RRL 
stars and whether it is associated to a metallicity gradient. 

To avoid spurious fluctuations when moving from the in¬ 
nermost to the outermost regions we ranked the entire sam¬ 
ple of RRL stars as a function of the radial distance. Then we 
calculated the running average by using a box starting with 
the first 200 more central objects in the list. We estimated 
the mean Rhb^, radial distance and standard deviation in 
each successive sample of 200, removing five objects in the 
ranked list and adding five, until we reached the last 200 ob¬ 
jects in the sample (651 RHBs, 536 RRLs, and 993 BHBs). 
We performed several tests changing both the number of ob¬ 
jects included in the box and the step size and found that 
our conclusions are not affected. The data plotted in the 
top panel of Fig. 5 show that the value of Rhb parameter is 
steadily increasing over a substantial fraction of the body of 
the galaxy (r<25 arcmin), but flattens to RHB ~0.55 in the 
outermost regions. The outermost regions display an almost 

^ We adopted the /j£fB=(B—R)/(B-|-V-|-R) parameter (Lee 
1990), i.e., the ratio among stars that are located either inside 
(V=RRLs), or to the blue (B=BHBs), or to the red (R=RHBs) 
of the RRL IS. 
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Figure 5. Top. The Rhb parameter as a function of the the 
elliptical radius. The vertical arrows mark the core (rc=5.8 
arcmin, Mateo 1998) and half-light radius (ri,=11.3 arcmin, 
Irwin &; Hatzidimitriou 1995). The vertical red lines display the 
errors on Rhb based on Monte-Carlo simulations. The Bottom. 
Same as the top, but for the ratio between bright and faint RRLs 
in the full sample. 

constant HB morphology {Rhb ~0.55). The steady increase 
in the Rhb suggests that we move from a HB mainly dom¬ 
inated by RRL and RHB stars {Rhb ~-0.1) inside the core 
radius (rc) to a HB mainly dominated by BHB stars and 
RRLs {Rhb ~0.4) beyond the half-light radius (rh). 

To further constrain the above intriguing evidence we 
also tested the radial trend of the ratio between Bt- and 
Ft-RRLs {clean sample), following the running average pro¬ 
cedure with a sample size of 50, replacing five objects in the 
ranked list until we reached the last 50 objects in the sample 
of RRLs. The data plotted in the bottom panel show that the 
above ratio is quite constant {NBt/Npt ~0.5) inside the rn, 
but it increases by a factor of three in the outermost regions 
(r~25 arcmin). We cannot reach firm conclusions regarding 
a flattening in the outskirts of the galaxy (r~30 arcmin) due 
to the limited sample of RRLs in these regions. The above 
findings therefore suggest that the metallicity gradient was 
already present >10 Gyr ago. 


6 CONCLUSIONS 

We have performed detailed multi-band homogeneous opti¬ 
cal photometry of the Sculptor dSph. We accumulated more 
than 4,400 images (proprietary plus archive images) cov¬ 
ering a time interval of 24 years. The precision of the 
photometry and the accuracy of both the absolute and rela¬ 
tive zero-points allowed us to provide new mean magnitudes 


for 574 variable stars. Among them 536 have been identified 
as Sculptor RRLs: 82 are new identification while for 219 
candidates we provide for the first time the pulsation pa¬ 
rameters. The new RRL data set allowed us to investigate 
several interesting open problems. 

i) True distance modulus. The use of three photometric 
bands and of two diagnostics, the PWRs in V, B~V, and 
in V, B-I, that are reddening free and only minimally af¬ 
fected by metal content, allowed us to provide very accurate 
distances. The errors on individual distance modulus are 
never larger than 0.02sys and O.OQran. To investigate pos¬ 
sible systematics affecting the quoted diagnostics, we have 
adopted three independent calibrations, named empirical, 
semi-empirical, and theoretical. We ended up with a mean 
distance moduli of 19.62 with cr^ = 0.04 using the clean 
(290) sample. All the estimates presented in this paper agree 
quite well with similar evaluations based on reliable distance 
indicators (Pietrzynski et al. 2008). 

ii) Metallicity range. We used the /-band PLR of RRL 
to constrain the metallicity range covered by Sculptor RRLs. 
The clean RRL subsample suggests that the Fe content of 
the old (t>10Gyr) stellar population ranges from -2.3 to 
-1.5 dex. This range is in excellent agreement with spec¬ 
troscopic measurements of RRLs, but signihcantly narrower 
than the distribution of RGB metallicities (-3.0 < [Fe/H] < - 
0.8 dex, Kirby et al. 2009). 

Hi) HB morphology. The use of homogeneous and accu¬ 
rate photometry covering a significant fraction of the body 
of the galaxy allowed us to constrain the HB parameter 
as a function of radius. We found that the HB morphol¬ 
ogy changes from -0.1 (dominated by BHB stars and RRL) 
to -1-0.55 (dominated by RHB stars and RRL) when mov¬ 
ing from the innermost to the outermost galactic regions, 
as expected from the presence of a younger and/or more 
metal-rich stellar population in the center of the galaxy (e.g. 
de Boer et al. 2012). 

iv) RRL luminosity distribution. We split the RRL sam¬ 
ple into Bt-RRL {V< 20.13 mag) and Ft-RRL {V> 20.13 
mag) sub-samples. The large sample of RRLs allowed us to 
constrain the radial gradient of the Bt/Ft RRL ratio. We 
found that this ratio is almost constant in the innermost 
regions while it shows a steady increase for radial distance 
ranging from rnto 25 arcmin. In passing we note that this is, 
to our knowledge, the first clear evidence in which HB stars 
and RRLs display similar radial trends. Analogous empiri¬ 
cal evidence was found in the Tucana dSph (Bernard et al. 
2008). However, a dichotomous population of RRLs was 
clearly identified there. 

This work shows that RRLs provide significant con¬ 
straints on the chemical enrichment history of its host 
galaxy. In particular, the combination of a large sample 
- one of the largest RRL samples in a dSph galaxy to 
date (Bersier & Wood 2002; Monelli et al. 2012) - with high 
quality, multi-band photometry allowed us to determine the 
metallicity distribution of the old stellar component, and 
revealed the presence of a significant metallicity gradient 
among the oldest stars (>10 Gyr). Our results are in agree¬ 
ment with the star formation history of Sculptor recently 
obtained by de Boer et al. (2012), the synthetic HB mod¬ 
elling of Salaris et al. (2013), as well as the various works 
based on spectroscopic data (see, e.g. Tolstoy et al. 2004; 
Battaglia et al. 2008; Kirby et al. 2009). In particular, we 
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find that Sculptor underwent a very efficient chemical en¬ 
richment: the metallicity range of RRLs shows an increase 
in the Fe content by almost one dex in less than 2-3 Gyr. The 
metallicity gradient suggests that star formation lasted until 
more recently in the inner regions, therefore further enrich¬ 
ing the interstellar medium. However, we cannot exclude the 
possibility of a merger at an early epoch, or tidally-induced 
heating of an initially disk-like structure into a spheroid 
(Mayer et al. 2001, 2006) as the origin of the gradient. 
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